Pointer States via Engineered Dissipation 
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Pointer states are long-lasting high-fidelity states in open quantum systems. We show how any 
pure state in a non-Markovian open quantum system can be made to behave as a pointer state by 
suitably engineering the coupling to the environment via open-loop periodic control. Engineered 
pointer states are constructed as approximate fixed points of the controlled open-system dynamics, 
in such a way that they are guaranteed to survive over a long time with a fidelity determined by 
the relative precision with which the dynamics is engineered. We provide quantitative minimum- 
fidelity bounds by identifying symmetry and ergodicity conditions that the decoherence-inducing 
perturbation must obey in the presence of control, and develop explicit pulse sequences for engi- 
neering any desired set of orthogonal states as pointer states. These general control protocols are 
validated through exact numerical simulations as well as semi-classical approximations in realistic 
single- and two- qubit dissipative systems. We also examine the role of control imperfections, and 
show that while pointer-state engineering protocols are highly robust in the presence of systematic 
pulse errors, the latter can also lead to unintended pointer-state generation in dynamical decoupling 
implementations, explaining the initial-state selectivity observed in recent experiments. 

PACS numbers: 03.67.Pp, 03.65.Yz, 07.05.Dz 



I. INTRODUCTION 

Understanding how and to what extent robust classi- 
cal properties can dynamically emerge in open quantum 
systems as a resuh of the unavoidable interaction with 
the surrounding environment has long been identified 
as a problem of central significance to quantum physics 
and quantum engineering. The concept of pointer states 
(PSs) has been introduced by Zurek [ij to capture the 
fact that not all initial pure states of an open quantum 
system may be equally fragile with respect to the interac- 
tion with the environment: PSs are distinguished by their 
ability to persist with high fidelity over time scales of 
practical interest, and are thus natural candidates to de- 
scribe "preferred" states in which open quantum systems 
are found in reality. As a result, PSs play a fundamental 
role in investigations of the quantum-to-classical transi- 
tion and quantum measurement models as well as of 
general aspects of "most classical" minimum- uncertainty 
states in quantum-dynamical systems 0-0]. In the con- 
text of quantum information processing, a set of mutually 
orthogonal PSs (a pointer basis) provides the simplest ex- 
ample of an "information-preserving structure" (IPS) 0: 
since arbitrary convex mixtures of PSs arc preserved, a 
pointer basis naturally realizes a robust classical memory. 
As a result, PSs are also practically attractive in view of 
their potential for long-lasting storage capabilities. 

From a physical standpoint, the robustness of PSs can 
be traced back to the fact that they become "least en- 
tangled" with the environment in the course of the dy- 
namics Mathematically (in a sense that will be 
made more precise later), this is only possible if the 
open-system Hamiltonian exhibits a sufficient degree of 
symmetry, which effectively allows PSs to be eigenstates 
(fixed points f^) in the resulting system-plus-bath di- 



lation. This has two implications: On the one hand, 
for a generic open quantum system, such a symmetry 
is not typical and at best approximate, thus PSs need 
not exist - with all the initial preparations of the sys- 
tem being rapidly degraded over comparable time scales. 
On the other hand, even in situations where a robust 
set of states would be naturally "cin-selected" over time, 
the latter would be inflexibly determined by the Hamil- 
tonian under consideration - with the resulting PSs not 
necessarily coinciding with states of interest, and with 
no control over their actual fidelity and lifetime. By re- 
versing this logic, we may then ask whether by suitably 
"engineering dissipation", that is, by resorting to external 
manipulation and perfection of the symmetry in the con- 
trolled open-system Hamiltonian, it is possible to make 
any target set of initial pure states into artificial PSs, 
so that the desired initial preparations can be robustly 
stored over time and their features retrieved on demand. 
This is the question that will be addressed and construc- 
tively answered in this paper - a task to which we refer 
to as "pointer state engineering". 

Our strategy relies on open-loop (feedback-free) quan- 
tum control methods, close in spirit to dynamical de- 
coupling (DD) approaches to decoherence siippression 
jsMlal and robust quantum computation |17H2G |. The 
idea is to start from the "bare" open-system Hamilto- 
nian that describes the interaction between the system 
and its environment, and to incorporate this Hamilto- 
nian along into a pre-dcsignated control recipe that acts 
directly only on the system, either in the form of suf- 
ficiently fast sequences of control pulses or continuous 
time-dependent modulation. Since in practice only lim- 
ited knowledge of the bath degrees of freedom may be 
available, an important requirement in DD constructions 
is that they be robust against variations in the bath op- 
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erators, and determined only by algebraic properties of 
the underlying open-system Hamiltonian. While a vari- 
ety of different protocols exist (and are being success- 
fully tested in the laboratory [2l], H^]), their common 
aim is the synthesis, perfection, and upkeep of an ef- 
fective Hamiltonian which is "dynamically symmetrized", 
so that interactions with the environment arc removed, 
up to a given order of accuracy, for the evolution of 
any state inside the entire system's Hilbert space [1, 
or within a (control-dependent) "dynamically generated" 
decoherence-free subspace or noiseless subsystem [23l[2^ . 
Technically, this is achieved via a number of algebraic 
techniques for manipulating generic interaction Hamilto- 
nians, along with analytic approaches for perturbatively 
(or numerically |25M27l | ) reducing the unwanted decoher- 
ence contributions and bounding the residual errors. 

While many of the building blocks used in the design 
and analysis of DD protocols will also be employed for the 
task of engineering PSs, it is important to clarify how the 
two problems differ. In a typical DD setting, the goal is to 
synthesize, for a given system-bath Hamiltonian, a target 
unitary propagator (the identity evolution) on the sys- 
tem with a sufficiently high gate fidelity, so that, ideally, 
arbitrary initial preparations can be robustly preserved 
over a desired storage time (or stroboscopically in time if 
the DD cycle is repeated). As a consequence, in a good 
DD scheme the control performance should be as unbi- 
ased as possible with respect to different initializations, 
and DD protocols must be constructed without making 
reference to and/or taking advantage of possible knowl- 
edge of the system's initial state. While the system-bath 
Hamiltonian is also given in the PS-engineering problem, 
the target set of initial preparations to be preserved is 
also specified as an additional input. Thus, knowledge 
of this target set should be explicitly incorporated into 
designing a good PS-protocol, so that the output fidelity 
is optimized for the desired PSs over the desired storage 
time. These differences result in important practical ad- 
vantages: while DD methods discovered so far can in fact 
only guarantee a reduced fidelity decay rate over time, in 
PS-engineering we will be able to guarantee a high fi- 
delity value over long time-spans. The end result is that 
the designated states can live much longer (in principle as 
long as desired) in a PS-cnginecring scheme rather than 
in general DD procedures intended for quantum memory. 

A few remarks may be useful to further place our 
work in context. First, we note that, interestingly, high- 
order/uncanceled corrections to effective Hamiltonians in 
DD constructions may contain an additional degree of 
symmetry (with respect to the minimum needed to en- 
sure the intended averaging of the systcm-bath interac- 
tion) and thus result in the generation of "accidental" 
PSs " and effective decoherence freezing, as predicted in 
d^ll^. Remarkably, such accidental PSs have been ob- 
served in recent DD experiments [22|,[33l ■ While the phys- 
ical mechanism responsible for the observed initial-state- 
scnsitivity is different than in a PS-cnginccring protocol 
(stemming, in the experiment, from systematic pulse er- 



rors), we shall also address the emergence of such acci- 
dental PSs in DD sequences and show how they may be 
formally related within a unified control-theoretic frame- 
work. Second, we also iterate that DD has been in- 
voked as a general strategy for generating IPSs in (non- 
Markovian) open quantum systems - decoherence-free 
subspaces and noiseless subsystem [l^, [l^, which can 
both be seen as multi-dimensional generalizations of PSs 
in an appropriate sense Q. Again, a key difference is 
that such schemes are not tailored to preserve a tar- 
get subspace or subsystem specified in advance. In this 
sense, our present analysis shares some motivation with 
(closed-loo p) s tabilization protocols for Markovian evo- 
lutions [3ll - [33| , restricted however to purely unitary con- 
trol resources and to the (simplest) case of discrete sets of 
states {classical IPS) in the system Hilbert space. Lastly, 
we recall that the possibility of a prolonged high-fidelity 
regime - so-called quantum fidelity freeze - has been ex- 
tensively analyzed in the context of Loschmidt echoes for 
closed quantum systems evolving under a vanishing time- 
averaged perturbation |34l436l and, in turn, shown to be 
intimately related to DD [33 ■ From this point of view, 
PS engineering may be interpreted as dynamically pro- 
viding the requirements for a stronger freeze phenomenon 
to emerge, in generic open quantum systems, for the tar- 
get set of initial preparations of interest. 



The content of this paper is organized as follows. In 
Sec. ini we describe the general control setting, along 
with the algebraic conditions required for the effective 
controlled Hamiltonian to admit PSs. In Sec. IIIII we 
characterize the effect that the unavoidable deviation of 
the actual effective Hamiltonian from the intended PS- 
supporting form implies on the quality of the engineered 
PSs. Our main result is a quantitative lower bound for 
the minimum fidelity which can be guaranteed for the 
PSs over a range of storage times in terms of both spec- 
tral/locality properties of the open-system Hamiltonian 
and details of the applied control scheme. Construc- 
tive protocols for engineering an arbitrary set of PSs are 
described in Sec. IIVI whereas explicit illustrations in 
paradigmatic control scenarios are given in Sec. IVl based 
on both semiclassical analytical results and exact numeri- 
cal simulations. In particular, single-qubit sequences are 
analyzed and contrasted to DD sequences resulting in 
accidental PSs, and sequences for engineering Bell states 
as PSs in two exchange-coupled qubits are presented, re- 
covering and extending partial results in |38| . Sec. IVII 
addresses the impact of limited control as resulting from 
imperfect initialization capabilities and/or imperfect con- 
trol operations. We both characterize the degree of ro- 
bustness of the PS-engineering sequences against system- 
atic control errors, and, conversely, show how such con- 
trol errors can result in accidental PSs in realistic DD 
sequences. Sec. I Villi concludes with a summary of our 
main results, along with a discussion of their possible 
implications and further open problems. 
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II. NOTATION AND PROBLEM SETUP 

Throughout this work, wc focus on a finitc-dimcnsional 
open quantum system, that is, a distinguished subsys- 
tem 5' with associated Ds-dimensional Hilbert space Hg, 
coupled to an environment (or bath) B, with associated 
Ds-dimensional Hilbert space T-Ib [H]. We use \j) to 
refer to system states and |6) to refer to bath states, the 
overall composite Hilbert space being given hy 'Hs^Si'Hb- 
For a product state \j) (g) \b), we refer to \j) (|6)) as the 
system (bath) component, respectively. A factorizablc 
density operator will have system and bath components 
in a similar way. Unless otherwise stated, we choose units 
in which h = 1. We also tend to drop trivial tensor prod- 
uct signs, identity operators on S {Is) or B (Ib), and 
ignore ordering of commuting operators as long as there 
is no ambiguity. Finally, we use the asymptotic notation 
with respect to operators in a lax manner. For example, 
A + 0{e) might refer to an operator A plus corrections 
whose norm (such as the maximum singular value norm) 
is 0(e). In perturbation theory, 0{x) will symbolically 
denote corrections that scale linearly with a quantity x 
which need not be dimcnsionless. 

The joint evolution of the system and the environment 
is taken to be generated by a Hamiltonian of the form 

Hit) = Ho + H,trm. 

Here, iJctri(i) is a controllable Hamiltonian that acts 
non-trivially only on "Hs and that shall be employed for 
generating unitary operations, whereas the free Hamilto- 
nian Hq is specified in a frame where no explicit time- 
dependence is present and can be further expressed as 

Ho = Hsb+Is®Hb (1) 
= Hs®Ib + H'sb + Is®Hb. 

That is, Hb denotes the internal evolution of the bath, 
and HsB denotes the interaction Hamiltonian between S 
and B, including the internal evolution Hs of S. For- 
mally, the latter can be isolated by demanding that H'g^ 
involves only traceless operators on i3. In this way, the 
limit of a closed (unitarily evolving) quantum system is 
recovered for H'g^ = (thus Hsb = Hs). The limit 
where the bath is treated semi-classically corresponds in- 
stead to Hb — C and the net effect of Hsb is a random 
modification of the system Hamiltonian, hence formally 
Hsb acts on Hs only [40| . 

Our goal is to use open-loop control, specifically 
through the application of a prc-detcrmincd Hctriit), to 
enhance and maintain the fidelity of an arbitrary target 
set of orthogonal PSs. The strategy we follow is to modify 
the evolution of the composite system by subjecting it to 
repeated identical control cycles. Let the evolution prop- 
agator for each period (cycle) of the controlled evolution 
be denoted by Uc- Our fundamental assumption is that 
the cycle propagator Uc, no matter how implemented, re- 
mains the same for all cycles applied. Suppose that the 



duration of each cycle is Tc and N repetitions are imple- 
mented, up to the total time T = NTc, with N integer. 
Then the cycle propagator Uc defines an effective cycle 
Hamiltonian He on Hs®Hb with the following structure: 

Uc{Tc) = Uc ^ cxp[-iTcHc], 

He = HsB,c + Is(® Hb.c, (2) 

where, in analogy with Eq. ([T]), Hb,c acts as an internal 
bath Hamiltonian and Hsb.c as a system-bath interac- 
tion, respectively. Wc emphasize that Hsb,c and Hb,c 
are not, in general, the same that appear in the bare in- 
teraction in Eq. ([1]). In fact, Hsb,c can be made much 
smaller than Hsb in DD, whereas pure-bath terms may 
be induced by the control even for a "non-dynamical" 
bath, for which the bare Hamiltonian Hb = 0. The 
periodicity of the evolution allows us to obtain a time- 
independent effective Hamiltonian He, since stroboscop- 
ically the evolution is given by 

Un{T) = Un^ {Ucf - exp[-zTi7j. 

The key to achieve our control objective is to design 
Hcf[\{t) in such a way that Hsb,c takes a special form 
that we define next. The actual control schemes that 
result in this special form will be given in Sec. IIVI 

Let the orthonormal set of states {|j)}^^Q [p < Ds) 
denote our target PSs. The basic step is to ensure that 
HsB,c is expressible in terms of a dominant operator 
Hdoni and a perturbation eifpcr, such that 

HsB.c{^) = -ffdom + ^Hpci-, 

ffdom = (3) 

where H' annihilates Span{|j)^^g} and Bj are distinct 
but otherwise arbitrary operators on B. This depen- 
dence of bath operators upon j is crucial to ensure that 
in the limit where e — )■ 0, each \ j) (but no coherent super- 
position) is invariant under evolution generated by He- 
Accordingly, we shall also refer to Eq. ^ as the "PS con- 
dition" henceforth. For nonvanishing e, the term ei/por is, 
in general, a system-bath Hermitian operator. Without 
loss of generality, we can shift any diagonal contribution 
of i/per in the PS basis into i?dom, so that 

{j\HpcAj)=0, Vj. (4) 

Note that if all but one PSs appear in the sum in Eq. 
([5]), the remaining state will automatically become part 
of the sum also, implying the equivalence of p = Ds — 1 
and p = Ds- If p < Ds — 1, the system Hilbert space 
is decomposed into two orthogonal subspaces, one gener- 
ated by the PS set {|j)}^^Q, and the rest. Beyond Eq. 
© and requiring e to be sufficiently small, we need not 
specify the controlled effective Hamiltonian further. 

Clearly, if the perturbation eifper was zero, and the 
system was initialized in any of the PSs, each of these 



4 



states would be indefinitely preserved with maximal fi- 
delity over time, just as the energy eigenstates in a closed- 
system setting. Realistically, however, even with sophis- 
ticated control schemes, the correction ei/pcr cannot in 
general be avoided. A good control scheme should en- 
sure that the states {|j)}j^o ^^^^^ singled out for their 
high-fidelity evolution, and can thus play the role of en- 
gineered PSs. Prior to providing explicit control schemes 
that synthesize effective Hamiltonians close to ffdom, we 
need to know how the corrections will affect the survival 
of the desired PSs in order to be able to maintain their 
performance. We thus proceed to analyze the extent to 
which the fate of the engineered PSs is modified in the 
presence of the inevitable correction terms. 



metric for quantifying the distance between ps{0) and 
j3s(N) is the survival probability (or input-output fidelity 
0), given by 

fN = Tr[ps{0)psiN)] = {0\ps{N)\0), 



which, since ps(0) is pure, is simply related to the cor- 
responding Uhlmann's fidelity via Jn ~ (/w)^- "^he 
initial state |0) will oscillate as a function of N with fre- 
quencies given by TcLOj{e). In the absence of degeneracy 
(for ojq), the amplitude of oscillation is determined by the 
difference between |0) and |0(e)) and is thus controlled by 
e through the ratio of the perturbation to the dominant 
term. More precisely, upon expanding |0) in the {|A:(e))} 
basis and applying , the fidelity loss reads as 



III. FIDELITY DYNAMICS OF ENGINEERED 
POINTER STATES 



I -In = 1- 



'^^^l(0|j(e)>P 



(5) 



As a basic motivating example, consider a two-level 
system, in which -ffdom is perturbed by a fixed intra-level 
coupling eHpar- If the energy levels are non-degenerate, 
the eigenstates of the unperturbed Hamiltonian i?dom can 
be seen as perturbations of the eigenstates of the per- 
turbed Hamiltonian. Similar to the familiar Rabi prob- 
lem, the discrepancy between the eigenstates is controlled 
by the ratio of the transition term to the energy gap and 
sets the amplitude of the oscillations that will ensue. For 
a degenerate system, the unperturbed eigenstates can 
end up being very far from the perturbed eigenstates 
and thus oscillate with a large amplitude [i^, p. 194]. 
The eigenstates of the original unperturbed Hamiltonian 
Hdom thus play the role of PSs which retain their high 
fidelity (modulo some oscillation) under perturbation. 

In an open-system setting, we shall use a similar idea 
with some modifications. The eigenstates of a perturbed 
system-bath Hamiltonian will still be close to the un- 
perturbed eigenstates as long as there is no degeneracy, 
but the actual difference between the eigenstates is ex- 
panded over the composite Hilbert space which can sam- 
ple over many basis states. The fidelity distance will thus 
need to reflect the larger size of the composite system- 
bath Hilbert space. This distance argument can be mod- 
ified using an ergodic argument to yield a tighter fidelity 
bound for the pointer basis elements. 



A. Semi-Classical Environment 

Consider first a semi-classical setting in which, for- 
mally, Hb,c ~ and Hsb,c is a system operator. Let 
{|j(^))}fe!fo~^ denote the (perturbed) normalized eigen- 



states of HsbA^) = Hdo 
ujjie). Thus, lim^^o \j{e)) 



eHpci, with eigenvalues 
\j) and lim^^o Ltjj(e) ~ ujj, 



where, without loss of generality, we may let a;o(0) = 0. 

Focus for simplicity on the preservation of a single PS, 
say |0). The system starts at ps(0) = |0)(0| and after N 
cycles evolves to ps{N) = ps{0)U~'^ . A convenient 



In the regime where e is sufficiently small, the required 
probability overlap (note that |(j|j(e))P is the so-called 
local density of states [13]) can he estimated using stan- 
dard 1st order non-degenerate perturbation theory. 



\jie)) = Z]l\ |j)+6E 



("|gpcrb) 



l("|gpor|i)r 



yielding 



l(ob-(£))l'«{,^y(o|Hp„|j)|Vc^|, j^o! 

Upon substituting the above expression in Eq. ([S]) and 
retaining terms up to order O(e^), one finds 



1 - /tv « ^4sin2[iVrcWj/2] 



l(0|i^pcr|j)|^ 



Under the perturbative assumption that |(0|iJper|j)/wj| 
is sufficiently small for all j 7^ 0, this finally results in the 
desired (uniform) fidelity lower bound: 



/jv > 1 - 4e^L'5 max 



l^porb-)P 



(6) 



As noted, if a degeneracy exists in -ffdom, the maximum 
fidelity loss cannot be bounded by a smooth function of 
e, since |0(e)) and |0) can be far apart regardless of e and 
at some point the initial state may oscillate too far and 
become completely lost. 

In the absence of degeneracy, Eq. ([6|) points to an 
elementary yet remarkable feature: no matter how long 
the time passed since preparation, the fidelity loss for an 
approximate eigenstate is small and bounded. This is a 
crucial property of PSs and one that we will strive to 
reproduce in the general quantum case. As it turns out, 
the fidelity dynamics of an engineered PS in the presence 
of a quantum environment will closely follow the above 



5 



perturbative derivation. Furthermore, the semiclassical 
approximation may be of independent interest in reahstic 
settings. For example, under appropriate physical and 
time-scale conditions, the contact hyperfine interaction of 
localized electronic spins with the surrounding nuclear- 
spin bath in semiconductors can be modeled in terms 
of a static but inhomogencous magnetic field as long as 
ensemble measurements arc considered [13, 13-111] ■ Our 
semiclassical analysis transcribes perfectly to this setting. 

Let us focus, in particular, on a two-level system, and 
let cTq, a = x,y,z, denote the corresponding Pauli ma- 
trices. In the semiclassical limit, let 

HsB = h - a = b^cr^ + byCfy + b^a^, 

where b = [b^^by^bz) is sampled from a distribution of 
random vectors. Clearly, Hsb itself does not result in 
a preserved pointer basis or a preferred direction unless 
the distribution of b is anisotropic. The effective Hamil- 
tonian, on the other hand, can be engineered as 

HsB,c = hz{h)az + e[h.^{h)(j^ + hy{h)ay\, (7) 

where /iz(b)(7z is the dominant term designed to pre- 
serve {|0) , |1)}. The small parameter e and the functional 
forms of /).a;(b), hyih), and /iz(b) will depend on the con- 
trol sequence used for producing Hsb,c, as well as the 
details of the probability distribution. The calculation is 
straightforward for simple control sequences. In our no- 
tation, cji = 2/iz(b) and |(l|Fper|0)|2 = h^{h)^ + hyihf . 
Our result for /at translates into 



N 



, [h4h)^ + hy jh)^] Sin' [NnKjh)] 
e , i»j 



as long as hz{h) ^ 0. Starting from /q = 1, the fidelity 
loss is thus at most e^/iz(b)~^[/i^(b)2 4- ft.y(b)^]/4. Eq. 
([5]) describes the fidelity dynamics for a single realization 
of the classical random field, whereas the actual open- 
system fidelity is the average of /jv over the distribution 
of b. We shall use this approach in Sec. IV A 11 to obtain 
approximate analytical results to be used in conjunction 
with our simulations of decoherence due to a quantum 
spin bath. Classical phase noise will also be relevant to 
the discussion of DD experiments on electron spins of P 
donors in Silicon in Sec. IVIII 



B. Quantum Environment 

In the quantum regime, Eqs. (H])-© describe a bipar- 
tite system, evolving under a Hamiltonian of the form 



He = (H. 



dom 



B,c) 



per- 



(9) 



While the dynamics generated by the unperturbed com- 
ponent Hjj would affect a generic initial state prepara- 
tion, it does not affect the fidelity of initial PSs. Con- 
sequently, in order to probe the fidelity dynamics of the 



PSs only, the basic idea is to isolate the "PS-preserving 
dynamics" by effecting a transformation to a suitable in- 
teraction frame. As it turns out, this transformation will 
play a crucial role in our derivation, as it will provides the 
basis for applying von Neumann's mean ergodic theorem 
(MET) [see Appendix H]. 

Let Ud = exp(— iTc-ffzj) denote the propagator for 
the pure-bath and dominant interaction dynamics. The 
single-cycle propagator can then be factored as follows: 



Uc = exp{-ieE)UD, 

OO ^ 

eE — eTcHpcv / — :[^TcH per. TcHd]^ 



(10) 



where £'[^+1 refers to corrections of second and higher 
order in TcHpcr, bm are the Bernoulli numbers, and 
L4, i3]„j = [■ ■ ■ [A, B], ■ ■ ■], B], with B appearing m times 
j49| . Note that due to our assumption in Eq. E is 
purely off-diagonal in the PS set, modulo the higher order 
]t;[2+] terms. The propagator for N cycles thus reads 



N 



N 



where the unitary operator Un represents the propagator 
in the toggling frame generated by {Ud)^ and can be 
approximated by invoking a Magnus expansion [50l | : 



Un 

[2+1 II 
N I 



exp(-iriAr), = e^^N + e'^^w^', 

N-l 

Y^iUDr'-EiUor, (11) 

0{N^Ef)^0{T'\\Hp,,f). 



The time-discretization implied by the periodic evolu- 
tion streamlines the application of von Neumann's MET: 
the sum in Eq. (jlip projects E onto the commutant of 
Ud and, as long as iJdom is non-degenerate, to that of 
-ffdom- To implement this exphcitly, let \j} (g) \b}j = \ j,b) 
and iOj^b denote the eigenstates and eigenvalues of Hd, 
respectively, and let us expand E in the operator basis 
induced by \j,b), that is: 



E 



E E 

jl ,bi j2M 



bi;j2-b2 I 



where Ej^ .j^ = {ji ,bi\E | , ^2) \ ji , &i) (j2 , 62 1 ■ Fol- 
lowing Appendix |3 we can partition fi^' in Eq. pT|) 

into diagonal components in the |j, 6) basis, namely fij^ 
(projected onto the commutant), and off-diagonal com- 
ponents in the \ j, b) basis, namely (projected outside 

the commutant). The diagonal components D^]^ vanish 

as a consequence of Eq. ([4]) , since fi^' is related to iJper 
via commutators with Hdom which preserve the commu- 
tant structure [5l| . The off-diagonal components $7]^ , on 
the other hand, do affect the PS fidelities, but instead 
of growing linearly with N (or time T), they are kept 
bounded by the MET. This is the fundamental feature of 
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PS preservation that guarantees a non-trivial long-time 
fidelity behavior. We can obtain an expression for the ef- 
fective off-diagonal terms after averaging under the MET 
[Eq. in Appendix E]: 



E ■ 



— ptJ- c\^3l .bi ^32:^2) 



where each term £-^1,61^2,62 rescaled by a factor that 
is controlled by the inverse of the energy difference 
Wjj.bj — ujj2^b2 s-nd can be bounded independently of N. 
In summary, we can approximate the total propagator 
up to time T as: 



Un = exp(-^e^!]V)(C/c)^ + 0{e^T^\\Hp^ 



where 



W^nW < 



, A = min |a;,- 

sin(r,A/2)' n^j2 



Notice that transitions within the same \j) sector do not 
appear at all, since E has no matrix elements within the 
same |j) sector, ultimately due to Eq. 

We now proceed to obtain an approximate upper 
bound for fidelity loss of the designated PSs, along with 
conditions for its applicability. For simplicity, let us as 
before focus on an initial PS preparation in |0), that is: 



p(0) = |0)(0|®pb(0), 



(12) 



where pb(0) is an arbitrary state on B. The joint state 
after N cycles becomes 

p(N) = cxpi^ienj^) [ |0)(0| (E> Pb{N)'\ cxpiienjj) 
+ 0{e^T^\\Hp,,f), 

where peiN) = U^pb{Q)U^^ . To focus even more on 
the fidelity evolution of the initial state, we can further 
partition il]^ into parts that couple to |0) and the rest: 

= n^"^ + VL'§''\ Up to corrections oi 0{e^\\VLj;j\\^)., we 
can then write 



p{N) = exp(-ier!l^>) [ |0)(0| ® pb{N) 

+ 0{e^r-\\H,,,,f) 

+ 0(e2T2||i/p,jVsin2(r,A/2)], 



ex-plied}^) 



(13) 



where we have bounded the corrections due to fac- 
toring cxp(— iefi]^) « exp(— lefl'^^) exp(— iefJ^^'^') by 
0(e^||il;^|p) estimated by the final asymptotic term. 



Thus, the fidelity loss is governed by efi^ 



|0) 



Sj#o,6i,b2('^o.''iy.''2 + ^-c-)- Then we may bound 



Let £^10) = 



in 



|0>| 
N I 



< 



em 



10)1 



sin(r,A|o)/2)^ 



^10) = I^J,f'2 

j^O.bi ,b2 



Physically, the "relevant gap" A|q) is the minimum energy 
difference between the |0) sector and any other sector in 



Hd. In principle, the latter can be refined if detailed 
knowledge of the energy levels of i?dom and Hg is avail- 
able. We also note that ||£:|o>|| < ||£^J1_~ Till^fpcrll- 
Using the fidelity bounds in Ref. [52|, |5j] [see also Eq. 
(|Bip ]. we can connect the operator norm of the effective 
Hamiltonian to fidelity loss. For small e, we have: 

l-fN< II + 0{e^T^Hp,,f) + 0{e^Hp,,\\y A') 

= — iiEi — y — (14) 

sin(r,A|o)/2) ^ ' 

+ 0(e2r2||i/p,j2) + 0{e^Hp,,\\yA^), 

which we can weaken and approximate to: 



fN>l 



perl 



sin(T,A|o)/2) 



(15) 



Notice that the fidelity bound in Eq. ([Ti)) . is similar in 
structure to the semi-classical case bound from perturba- 
tion theory in Eq. ^ , but is a weaker bound by a power 
of 2. While this is due to the coarse bounds for fidelity 
loss in terms of effective Hamiltonian from Ref. [11], we 
emphasize that these bounds have the advantage that 
they apply irrespective of the state of the environment. 
Also notice that the appearance of operator bounds guar- 
antees a polynomial dependence on the number of subsys- 
tems ub in the environment, as all system-environment 
operators can be written as sums of few-body interaction 
terms among the environment subsystems. 

The bound on the fidelity loss in Eq. (|14p is valid 
provided that two conditions are met: first. 



eTWHr, 



« 1, 



(16) 



which is a technical requirement for our approximations 



to be valid (in particular, for neglecting fijy 

spect to rij^'), and delineates the regime of a] 
of MET in our analysis; and second, 



[2+] 



with re- 



A > e||iJ, 



perl 



(17) 



which is needed to ensure that no transition between 
the different pointer sectors occurs due to the pertur- 
bation. Both these requirements can be met by reducing 
e, the perturbation strength. Interestingly, we may inter- 
pret the gap condition in Eq. (jl7p as a simple energetic 
protection of coherence: PSs are protected indefinitely if 
the perturbation is too weak to cause a major transition 
among them. The role of a sufficiently gapped i?dom is 
precisely to ensure this energy barrier. 



IV. CONSTRUCTIVE CONTROL PROTOCOLS 

The results of the previous section provide a quantita- 
tive estimate on the maximal fidelity loss that PSs may 
incur, no matter how the latter happen to be produced. 
Ideally, the system starts with a nearly-perfect prepara- 
tion (by means of a suitable quantum operation) in one 



7 



of the PSs, thus the fidchty /g = 1. As the control cycles 
are repeatedly applied, the fidelity will decay quickly for 
an initial transient [see Appendix [B] , but will eventually 
enter a "saturation" regime with a fidelity level lower- 
bounded by Eq. (fTS)) once the number of applied cycles 
N is sufficiently large for the MET-averaging to kick in: 

NTcA = TA > 1. 

This saturation regime will hold as long as the conditions 
in Eqs. P^ - ((T7)) arc satisfied. Within this picture, to 
engineer an arbitrary set of PSs, we need to: 

• PSl: Engineer the dominant cycle Hamiltonian iJdom 
into the diagonal form of Eq. ([3]), so that the PS 
condition holds for the desired choice of {|.7)}^^0' ''^ith 
l<p<Ds-l; 

• PS2: Minimize the perturbation iJper, by enforcing 
high modulation rates (that is, by reducing Tc) and/or 
employing advanced control protocols similar to those in 
high-order DD [ll, [l^ 113] to get a higher fidelity bound 
[See Sec. IV C| . Notice that this also implies that Eqs. 
([T5|) and ([T7|) can be better satisfied and thus the desired 
PS set preserved to a longer guaranteed time. 

In this Section, we describe how the tasks PSl and 
PS2 can be achieved using unitary control pulses, for 
an arbitrary open-system setting in which Hsb and Hb 
are quantitatively unspecified (bounded) operators on 
Hs ® "Hb- Although in principle continuous modula- 
tion schemes could be envisioned for this purpose, the 
pulsed control setting we focus on has the advantage to 
be mathematically straightforward while providing inter- 
esting connections to both DD theory and recent experi- 
ments. We begin by assuming that Hsb is generic - that 
is, it has no special symmetry hence no degeneracies - so 
that the engineered PSs and only those remain at high 
fidelity. We revisit this assumption in Sec. IIVBI 



A. Generation of the Effective Hamiltonian 

The control cycles in our construction are pulse se- 
quences in which free evolution intervals of duration 
are punctuated by application of system unitary oper- 
ators Pi, so that Tc = The operators Pi are 
assumed to be implemented as (nearly) instantaneous 
pulses by Hcu\{t). While in this section as well as in 
the examples of Sec. |V] we will also assume each pulse 
to be implemented perfectly, it turns out that such ideal- 
pulse assumptions are not essential as long as the unitary 
propagator does not change from cycle to cycle [See Sec. 
IVI| . Similar to DD, the control pulses are so designed 
that Pi cancel each other: P„ • • • Pi = Is- The overall 
unitary propagator for the cycle is an ordered product 
of pulse unitaries interlaced with evolution propagators 
exp(— iTji/o) associated with the free intervals in which 



Ha = Hsb + Hb- That is, 

Uc{Tc) = P„exp(-^T„i^o)•••PlCxp(-^Tli/o) (18) 
= QI^ cxp(-iT„i?o)Qn ■■ - QI cxp{-iTiHo)Qi 
= cxp{-iTnQiHaQn) ■ ■ ■ cxp{-iTiQ\HoQi) 
= exp{-iTcHc), 

where Qj = Pj-i • • • Pi for n > j > 1 and Qi = Is- The 
effective cycle Hamiltonian He [Eq. ([2])] can be approxi- 
mated using a Magnus expansion: 

n 

TcHc = rjQ]HoQj + 0{T^cHl). (19) 

i=i 

In the special case of a uniform sequence, for which = 
T = Tc/n, we will simply denote the cycle in Eq. (jl8p by 

/Pl/P2---/P„, 

where operations are now applied left-to-right and / 
stands for a free evolution interval. 

The task PSl then reduces to casting He in the form of 
Eqs. (H])-® for a designated PS set {|j)}^Io. Note that 
in non-selective (universal) DD schemes @, a fixed choice 
of Pi (determined by the system dimensionality Ds and 
the symmetry of Hsb) is used to cancel iJo (modulo pure- 
bath terms) up to various orders in |(rc^^o|| ■ br this sense, 
PSl is easier than DD, as it requires cancellation of far 
fewer terms in H^. In fact, the minimum number of in- 
tervals in a cycle designed to completely cancel a generic 
Hsb is D'g 0, [5^, while in contrast, as we shall prove 
below, the cycle for achieving PSl for a complete pointer 
basis requires only Ds time slots. Achieving PSl (and 
PS2) is. rather, close in spirit to selective DD, the price 
to be paid, however, being that the required control op- 
erations become necessarily state- dependent. A similar 
idea has been explored in Refs. [ssl. |54|. [56j . where "ro- 
tated" DD pulses are used to "lock" states in a particular 
one- or two- dimensional subspace of a two-qubit dissi- 
pative system. Note that selective (and encoded) DD for 
subspaces [13, HI] has also been specifically invoked for 
dynamical quantum error suppression, in alternative to 
non-unitary control strategies based on error-correcting 
codes [5§| . A basic difference between the PS-engineering 
protocols that we will provide and the above-mentioned 
methods is the dimension of the preserved structures. 
Here, the preserved states correspond to isolated points 
(zero-dimensional manifolds) in the system Hilbert space, 
while by decoupling subspaces (or subsystems), the cor- 
responding IPS corresponds to continuous manifolds of 
higher dimension. 

As mentioned, the desired set of PSs may comprise any 
number p of orthonormal states. For clarity, we start 
with the task of engineering a single PS, follow with pre- 
serving a complete pointer basis, and finally provide a 
general recipe for engineering an incomplete pointer set 
that includes the former protocols as special cases. 
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1. Single Pointer State 

Let |0) be the solitary target PS {p — 1). Define the 
unitary reflection operator 

g = 2|0)(0|-/5. (20) 

Notice that despite the simple representation, implement- 
ing Q need not be simple, and we further comment on 
that in Appendix [Cl Consider a uniform pulse sequence 
consisting of two Q pulses separated by equal intervals of 
duration r = Tc/2. Using Eq. (fT9|) . one can verify that 
He for this sequence obeys the PS-condition in Eq. ([3]), 

= |0) (0| ® Bo + ^^pcr + Is®Hb, 

where Bq is a bath operator and iJpor — 0(r||_ffo|P) 
(ssl [s^ l- Since the latter can be made smaller by us- 
ing shorter pulse intervals, the sequence will preserve |0) 
with an arbitrarily high fidelity that can be maintained 
for a long time. Thus, PSl and PS2 are both achieved 
by the control cycle fQfQ = QQ- 

If Ds = 2, the choice of |0) determines a unique or- 
thonormal state |1) and Eq. ^ yields Q = |0) (0| - 
|1) (1| = (Tz, resulting in a sequence that we shall simply 
denote ZZ henceforth. In a DD problem, this sequence is 
designed to suppress terms of the form ax<E)Bx + o-y<E}BY 
in HsB- In contrast, a universal DD sequence (such 
as XYXY) is designed to cancel every possible term in 
Hs B , including cr^ Bz ■ This generally results in an en- 
hanced quantum memory where all states are preserved 
with a higher fidelity initially, but need not lead to a fi- 
delity which is maintained for any particular state: even 
approximate eigenstatcs induced by XYXY are not ap- 
parent. Despite this distinction, in Sec. IV Al we will ana- 
lyze the structure of Hsb,c for both sequences in detail, 
and show how the states |0) and |1) form a (accidental) 
pointer basis under the XYXY sequence as well. 



2. Complete Pointer Basis 

Let {|j)}j^o~^ denote the set of p = Ds states we de- 
sire to preserve. Following Issj . let us define the unitary 
operator as 

Ds-l 

aos= J2 ^'^'b'XjI. a. = e2-^/^^ (21) 
i=o 

Direct calculation shows that ioi < i, j < Ds — 1, 

Ds-l 

J2 i'^Dsy'^Mjl ® B,j){aDsf = DsS;Ji){i\ ® Bu. 

k=0 

Now consider a control cycle consisting of equal free in- 
tervals of length r = Tc/Ds, during which is ap- 
plied Ds times, or /(Tds ' ' ' I^Ds- Thus, Pj = and 



Qj = i'^DsV ■ Using Eq. ([T5|) . the effective cycle propa- 
gator again takes the desired form of Eq. : 

Ds-l 

Hc--Yl \j){j\®Bj+Hp^, + Is®HB, 

3=0 

where, as before, ||-ffpor|| = 0(t||7Jo|P)- 

Note that the above sequence is by no means uniquely 
suited for the tasks PSl and PS2. For example, sequences 
based on Hadamard arrays [63, HH or products of CTz 
Pauli matrices [g^] can also be envisioned when the sys- 
tem is a collection of Ug qubits {Ds = 2"''). Nonetheless, 
the sequence of aos operators described here is notable 
since it requires a single pulse type. As we shall see in Sec. 
IVIBl this implies an extra degree of robustness against 
systematic control imperfections. 

3. Incomplete Pointer Basis 

Let {|j)}^^g denote the set of p < Ds — 1 orthonormal 
states we wish to preserve (the case p = Ds — 1 is equiva- 
lent top ~ Ds, as already noted. The set {\j}}^Zo j possi- 
bly along with other orthonormal states {\j)}fjjZ^ , forms 
an orthonormal pointer basis for "Hs- The construction 
described in the previous subsection can be readily mod- 
ified to this scenario by introducing the following unitary 
pulse operator: 

-?)s - E-^'+'b->01 + \m, ^ = e2"/(^+i). (22) 

The 1st sum in Eq. corresponds to a diagonal block 
of size p while the 2nd sum implies that cr^^ acts as iden- 
tity on Span{|j)}j^p~^. Consider a control cycle consist- 
ing of equal free intervals of length r = Td {p+1), during 
which (7^ is applied p + 1 times, or /cr^ ■ ■ • /c^ . • Us- 
ing Eq. m with Q, = {a^y for j ^ 1,...,J>+ 1, 
we can verify that the effective cycle propagator for this 
sequence obeys the PS-condition in Eq. ([3]): 

p-i 

t/e = exp[-ir,(^ ® Bj +H' + Hp,, + Is® Hb)], 

3=0 

where, as required, H' annihilates Span{|j)}^^Q and 
again ||i7pej = 0{T\\H^f). 

B. Over-Symmetric Systems 

In essence, PS engineering is about introducing and 
maintaining symmetries in the open-system Hamiltonian 
Hq, and indeed the protocols presented in Sec. IIV Al in 
conjunction with the error bounds of Sec. IIIIl succeed at 
achieving this goal if Hsb is generic, that is, not over- 
symmetric. The occurrence of unwanted symmetries can 
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cause two distinct problems. First and most importantly, 
they can result in energy degeneracies in 77dom- This will 
affect the perturbation theory requirements in our deriva- 
tion of fidelity of PSs, since vanishing energy differences 
[ujj and ujj^b in Eqs. ([6]) and ([15]), respectively] will lead 
to a divergence in fidelity loss. In such cases, obliviously 
applying the above control procedures may result in PSs 
that decay quickly. Second, another problem arises when 
a state other than the designated set {|i)}^^Q happens 
to satisfy the PS condition in Eq. ^ due to additional 
symmetry. The intended PSs will then be preserved, but 
not exclusively so. In particular ^decoherence-free sub- 
spaces or more general IPSs [sl, Q could exist/emerge 
in the presence of the control sequence that is used for 
engineering PSs. 

Both the problematic scenarios of energy degeneracy 
and non-exclusivity of PSs can be remedied by modifying 
the control protocol so that Hq (and consequently He) is 
suitably "de-symmetrized". In practice, a situation of un- 
conditional failure of the protocols described in Sec. IIV Al 
hints at hidden symmetries in the bare Hsb that persist 
in the effective cycle Hamiltonian Hsb, a and need to be 
addressed by modifying the control protocol so that the 
transformation from Hq to H^ not only introduces the 
desired symmetry but also removes the undesired ones. 
A straightforward way in which such a desymmetrization 
can be obtained is to "turn on" an constant Hamiltonian 
on the system. For degenerate cases, an explicitly di- 
agonal and non-degenerate control Hamiltonian can be 
applied, whereas for the undesired PSs an explicit cou- 
pling to yet another state (if available) will guarantee the 
desired desymmetrization. Another (pulsed) solution is 
to adjust the previous control protocols to remove the 
over-symmetry, as we describe next. 

For simplicity, let us focus on a simple case where a 
single state |s) is the source of the problem: it either 
shares energy eigenvalues with another fellow orthogonal 
PS, or it is an undesired PS. We will assume that another 
state |i), orthogonal to \s) and the rest of the PSs, exists 
and that the following unitary (self-inverse) operator, 

R^\s){t\ + \t){s\+ J2 

is available for control. Notice that R permutes the states 
\s) and \t). Let the control cycle for generating the de- 
sired PSs [Sec. |IV] be given by fP ■ ■ ■ fP, where each 
pulse is originally applied n times. We will ammend this 
sequence by concatenating [l^ it with fRfR, that is, 

IfRfR] P[fRfR] P---[fRfR] P 

As a result, the number of pulses is now 2n and they 
alternate between RP and R. Using Eq. p^ . the cycle 
Hamiltonian becomes: 

HsB,c = P~\RHsbR + Hsb)P\ 

fe=0 



in which the inner expression (RHsbR+Hsb) effectively 
removes the degeneracy by "mixing" the states |s) and 
\t) as long as there is no further symmetry among them. 
Notice that, if the symmetry of Hsb is manifest, such as 
in a system exposed to collective noise, the operator R 
can be chosen from the outset to be a product of unequal 
unitary operations on each subsystem, which inevitably 
results in desymmetrization. 

V. ILLUSTRATIVE APPLICATIONS 

A. Single-Qubit Sequences 

We begin by considering the task of engineering two 
arbitrary orthogonal states |0) and |1) as a pointer basis 
in a qubit. This basis defines the relevant Pauli opera- 
tors, and the bare system-bath interaction Hsb can be 
correspondingly expanded as 

Hsb ^ (Jx <^ + ay (E) By + (g) B^, 

for generic (traceless) operators {Ba}- As before, let Hb 
be the bath internal Hamiltonian. Our goal is to thus 
preserve the eigenstates |0) and |1) of az and, as men- 
tioned, we shall both analyze the ZZ PS-protocol and 
the XYXY universal DD sequence. We can approximate 
the effective cycle Hamiltonian He using the Magnus ex- 
pansion. For the ZZ sequence we have Tc = 2r and 

hIZ^ = ^^®Bz, (23) 
eijiffi = U^ax®Bx + (Ty®By,(Jz®Bz + HB] 

+ 0(t2). (24) 

The perturbative parameter in this setting is clearly pro- 
portional to tIISqII, where a = x,y. For the XYXY 
sequence we have Tc — 4t and 

^d™""' = ^T(2tax®[Bx,HB] 

+ az®{i[Bz,HB]+{By,Bx})y (25) 

whilst ||i?pCT^'^'*^' II = O(r^). For both sequences, higher- 
order corrections in effpor will gcnerically contain pure- 
bath terms, causing Hb,c to differ from the bare Hb- 

Clearly. H^^^^ has separable eigenstates of the form 
\i) eg) \bi) for i = 0, 1. For hI^^^^\ these eigenstates 
appear independently of Bx only if Hb = 0, that is, the 
bath is non-dynamical. Physically, the non-dynamical 
bath is a special but important case which provides, in 
particular, a prevalent approximation for decoherence of 
localized electronic spins in semi-conductors [2^ [2^ . In- 
terestingly however, even for a generic dynamical bath, 
|0) and |1) are still preserved under the XYXY sequence. 
To see this we have to redefine -ffdom- Recall that the fi- 
delity preservation of PSs relies on the averaging caused 
by Ho, which since -ffdom is made arbitrary small by 
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shrinking r. is dominated by Hb- Consider the bath- 
operators appearing in H^^^ in Eq. (|25p . The key 
observation is to reahze that for any operator of the form 
[Ba , Hb] , the component afong the commutant of Hb 
vanishes (as long as Ba is traceless), and hence these 
terms will act solely as perturbations according to Eq. 

In contrast, the anti-commutator {By,Bx} will typ- 
ically have a significant component along the commutant 
of Hb, and will break the symmetry in favor of the cr^ 
eigenstates. Thus, we may redefine the dominant Hamil- 
tonian responsible for PSs as 



H 



IXYXY],\\ _ 



4Ta,(E){By,B^}\ 



(26) 



which clearly still preserves |0) and |1). A quantitative 
analysis of the resulting long-time fidelity follows. 



1. Semi- Classical Environment: Analytical Results 



The operators h]^^'^^^ 



and ^^dom '^^^ calculated 
analytically within the semi-classical approximation of a 
static random-field environment described in Sec. IIII Al 
Phenomenologically, such a random magnetic field b = 
{hx,by,bz) can approximate the effect of a nuclear spin 
environment interacting with a central spin qubit system 
under ensemble measurements (4^ . In addition, as we 
shall describe in Sec. I VII systematic pulse imperfections 
can also formally mimick a classical always-on magnetic 
field. By invoking Eq. (|5]) with h^ih) = bz, we obtain 
the ensemble-averaged fidelity loss: 



[hx{h)^ + hy{h)^] sin^ (NTrbz) 



(1 - /iv)b = ^^^P(b)db, 

(27) 

where P{h)dh is the probability density associated with 
the distribution of b. In the limit of od, the rapidly 
oscillating term sin^ [NTcbz) is smoothed out into a sat- 
urating (plateau) behavior. In contrast, for states other 
than the Cz-eigenstates, the error grows to a maximal 
(unit) fidelity loss value. 

It is illustrative to compute the average fidelity loss 
for the ZZ sequence with a specific distribution for b. 
The calculation carries over to other sequences and can 
easily incorporate models of pulse imperfection reflected 
in different distributions for b [see Sec. IVI| . Let us as- 
sume, specifically, that the distribution is isotropic and 
the magnitude b is distributed according to a normal dis- 
tribution P{b) with zero mean and standard deviation B. 
The expressions for hxih) and hy{h) can be read off Eq. 
(P^ . The average fidelity loss can thus be calculated in 
a straightforward manner. For small N , we have: 



(1 - /jv)b 



B'^N'^T^ + O (N'' 



1 

10' 



-B^T\^. 



(28) 



In this regime, the fidelity decays quadratically with both 
the elapsed time T and the control time scale t, also in 
line with general error bounds for cyclic DD [65|- In 



contrast, for large A^, the average fidelity loss saturates 
to a limiting (time-independent) value: 



2^2 



(1 - /oo)b = 



B't 
12 



0{e- 



2^2 



A^- 



(29) 



which is controlled by the product Bt. The pattern of 
initial fidelity decay, followed by saturation, is the same 
pattern observed in exact numerical simulations of DD in 
quantum dots [28l . [29[ . In such a setting, the bath consists 
of a large number of nuclear spins, each coupled to 
the central spin with a strength j„i a-nd with 



1/2 



riB 



being a measure of the coupling strength [see also the 
upcoming Eq. (|31]) ]. In the semi-classical limit, the en- 
vironmental spins produce a classical Overhauser field. 
Using this analogy and the exact result for the asymp- 
totic fidelity derived in Ref. [2^ [2^ for a maximally 
mixed bath initial state, that is, 1 — f^o = A^ub/^Q, 
we can interpret the variance of the i?-ficld as follows: 



B2 = -A^UB. 



(30) 



Thus, we recover the polynomial dependence of the long- 
time fidelity of the engineered qubit PS's on ub, as ar- 
gued in Sec. IIII Bl on general grounds. This will be also 
verified in the upcoming numerical results [Figs. [SHI]. 

The initial decay of fidelity followed by a saturation 
(or freeze) regime is consistent with the more general 
semi-classical analysis carried out in [33 - l36j . The semi- 
classical approach also highlights an interesting connec- 
tion to the Krylov-Bogoliubov method of averaging in 
classical dynamics [66j . where an oscillatory motion is 
replaced by an average that yields approximate integrals 
of motion. In the quantum treatment, a similar behav- 
ior emerges from the theoretical model developed in Sec. 
IIII Bl whereby the initial decay is followed by a bounded 
fidelity loss in the MET regime. Note that unlike the 
classical limit, after many oscillations (possibly when T 
no longer satisfies Eqs. or p7)) ). fidelity may even- 
tually drop, as the linear approximation implicit in the 
MET theorem (ignoring ilt'^+l terms) need no longer hold. 



2. Quantum Environment: Exact Numerical Results 

We illustrate our findings by means of exact numeri- 
cal simulations where a spin qubit [S = {Sx, Sy,Sz) spin 
operators] couples through Heisenberg interaction terms. 



HsB = ^ jmS • I' 



(31) 



m=l 



to a quantum environment consisting of spin-1/2 parti- 
cles [I^™) spin operators]. The (bare) internal Hamilto- 
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Figure 1: (Color online) Fidelity error (1 — /) as a function of time for XYXY and ZZ sequences as well as the free evolution 
for \+X) — (|0) + \ l))/y/2 and \+Z) = |0) initial states. Left: ns = 15 bath spins. Right: ub = 8 bath spins. In both cases the 
bath is non-dynamical (/3 = 0), and the pulse interval is set at r = 0.01 in arbitrary units. In order to allow for a comparable 
effective _B-field [Eq. (|30p ] and error range, the maximum coupling \ Jl\ ~ 1 {\Jr\ = 4) in the left (right) panel, respectively, 
with time kept in units of 1/ Jl in both cases. The difference in the time spans is due to practical limitations. Note that here 
and in the figures that will follow, the time axis is sampled non-uniformly to emphasize the interesting features of each time 
series without using too many data points. The sample points are connected by lines as a visual aid. 



nian of the spin bath is taken to be of dipolar form: 

riB 
rn— 1 k<m 

The coupling constants jm between the system spin and 
each bath spin are arbitrarily generated by (uniformly) 
randomly sampling between — J and J. Similarly, the 
coupling constants Pmk between each two bath spins are 
randomly chosen between — /3 and /3. We also assume 
that initially the bath is in a fully mixed state. 

Fig. [1] compares the fidelity loss of different initial 
preparations under free evolution and under the ZZ and 
XYXY sequences designated to preserve |0) = \+Z) and 
|1) = \—Z). Specifically, the system is prepared either as 
\+Z) (intended as a PS) or \+X) = (|0) + |l))/\/2 {not 
intended as a PS). We observe that for the free evolution 
or for state |+^), the fidelity loss approaches a maximal 
value corresponding to a completely mixed state. Degra- 
dation is significantly faster for the free evolution and 
ZZ with l+v^C)) and considerably slower for the XYXY 
(with l+X)), as DD will nonetheless universally extend 
coherence times. In fact, the state \+X) in the model 
with the larger environment (left panel) did not reach 
the maximal loss state within the simulation time under 
the XYXY sequence. On the other hand, with either 



the ZZ or XYXY sequence, the desired PS \+Z) is pre- 
served with high fidelity after the initial transient is over. 
This initial decay is expected to occur before the MET 
regime kicks in and is approximately described by Eq. 
(PS)) . Notice that the value of fidelity is non-trivially 
bounded after the initial decay, and furthermore it con- 
verges smoothly to a saturation value for all engineered 
PSs in Fig. [T] The saturation behavior is in agreement 
with the semi-classical result in Eq. ([2^ . We also point 
out that the initial fidelity loss values depicted in Fig. [T] 
correspond to a the first period of control, at T = Tc, 
which is expected to scale with ||TciJper|p. 

Fig. [5] depicts the change in fidelity of the preserved 
state |0) by applying the sequences at different pulse in- 
tervals r. As expected, faster modulations (smaller r) re- 
sults in a smaller perturbation strength and thus higher 
fidelities. For both sequences, the ratio ||ei?por/-ffdom|| = 
0(t) determines the long-time fidelity loss [Eq. (fTS)) ]. but 
the difference in the structure of iJpor vs. H^om in the 
two cases [Eqs. ([^ - (^5)) ] is responsible for the notice- 
able difference between the final fidelity saturation value, 
the ZZ sequence outperforming the XYXY DD protocol. 

In Fig. [3l we directly probe the dependence of the 
limiting fidelity (1 — /oo) on the pulse interval r of the 
applied sequences, for various randomly generated cou- 
pling patterns {jm} [Eq. (pij) ]. We can readily verify 
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Figure 3: (Color online) Long time (saturated) fidelity error 
(1 — /oo) as a function of rescaled interval tA for preserving 
|0) with the ZZ sequence, and A = (Em'Li Jm/^s)^''^- The 
expression 1 - /oo = 0.95(ry4)^-33 [not shown in the fi gure] is 
a power-law fit to the series with the smallest A. The number 
of bath spins ub = 8 and Hb = 0. The overlapping curves 
correspond to randomly generated coupling patterns between 
the system qubit and the bath spins. 



Figure 2: (Color online) Fidelity error (1 — /) of ]0) as a 
function of time for XYXY and ZZ sequences with different 
pulse intervals r. The number of bath spins ub = 15 and the 
coupling is capped at | J| = 4. 



that the observed fidelity loss scales with (At)'^, as ex- 
pected. In addition, Fig. |4] superposes the fidelity satu- 
ration values corresponding to environments of different 
size {ub = 3, • • • ,8), with the horizontal axis set to the 
bath-size- rescaled pulse interval Arn^^^ The fact 

that there is little variation in the curves confirms the 
validity of the semi-classical approximations leading to 
in Eqs. (| 2 9 1) - ((30)) . and provides evidence that our sim- 
ulation results may be reliably extrapolated to realistic 
environments with a far larger number of spins. 

The numerical data presented thus far have addressed 
only "non-dynamical baths", in which the internal cou- 
pling strength parameter /3 = 0. In the more general 
case where /? 7^ 0, the fidelity saturation value shows no 
significant dependence on /3 over the parameter range we 
explored. Results are summarized in Fig. [S] The inde- 
pendence of saturation fidelity from /? is expected within 
the validity of the MET regime as well as in the semi- 
classical approximation. On the other hand, the value of 
fidelity loss right after the initial cycle does depend on /3, 
and more strongly so for XYXY than ZZ. 

Finally, note that for initial states other than \X+) 
that are closer to \ZzL) PSs, the fidelity loss is smaller 
than the fidelity loss for In Fig. [SJ we compare 

the relative survival of initial states in the xy-plane of 



the Bloch sphere, IV') = cos(6'/2) |0) -h sin(6l/2) |1), for 
various values of 9 for free evolution and the sequences 
XYXY and ZZ, respectively. We have defined the rel- 
ative survival of a state as the inverse ratio of its long 
time fidelity loss and the maximal fidelity loss, that is, 
am = (1 - /min)/(l - /oo(l^))), where 1/2 for 

a qubit (corresponding to the fully mixed state) . Clearly, 
the engineered PSs {9 = 0) enjoy the highest survival ra- 
tios for sufficiently small r's [Eq. (|29)) ]. as intended. Also 
notice that significant difference between the engineered 
PSs and the other states is virtually non-existent in the 
free evolution. 



B. Bell-State Pointer Engineering 

We next consider a two-qubit spin system, coupled to a 
spin- 1/2 environment similar to the one described in the 
single-qubit case. The spin operators for the two qubits 
are now denoted by S^^^ and S^'^\ respectively. Each 
qubit interacts individually with the bath spin particles 
via a Heisenberg Hamiltonian, thus 

k—1 m— 1 

where as before we take the jm,k coupling con- 
stants to be sampled from a random distribution with 
maxm.fc lim.fcl < J, and no additional symmetry is 
present (in particular, ^ jm.2 for at least one m). 
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Figure 4: (Color online) Long time (saturated) fidelity error 
(1 — /oo) as a function of bath-size-rescaled interval rAn^J'^ 
for preserving |0) with the ZZ sequence. The simulation used 
ub = 3, • • • , 8, with randomly generated couplings. The re- 
maining parameters are the same as in Fig. [S] 



The bath internal dynamics is governed by the dipolar in- 
teraction described by Eq. ([5^ . where we choose /3 = J. 
In addition, we also allow for an always-on Heisenberg 
interaction to be present among the qubits, that is, 

For concreteness, unless otherwise specified we shall 
match the intra-qubit interaction strength to that be- 
tween the qubits and the environment spins, K — J. This 
corresponds to a coherent evolution between the qubits 
that is comparable to (and competing with) the interac- 
tion with the nuclear spin environment. We will use pulse 
sequences described in Sec. II VI for preserving various sets 
of the Bell basis states: 

= (|01) + |10))/V2, 
= (101) - \W))/V2, 
= (|00) + |11))/V2, 
= (|00)-|11))/V2. 

In particular, we engineer three sets of PSs with p = 
1,2,4 states chosen from the Bell basis. In what follows, 
all operators are given in the computational basis. 

The sequence El is designed to preserve |EPRi) (sin- 
glet state) only. The cycle for El is based on the pre- 
scription given in Eq. pO)). with |0) being chosen as the 
EPR basis elements. With respect to the computational 
basis, the required control cycle consists then of two ap- 
plications of the SWAP gate (see also [38|): 




Uei = 



/ 1 
10 
10 

Vo 1 



r/4 




- o - 


XYXY, = 


= 


- e - 


XYXY. I'i = 


= 0.1 


- o - 


XYXY. I'i = 


= 1 


- - 


XYXY, 13 = 


= 1(1 




ZZZZ, 13 = 







ZZZZ, 13 = 


0.1 




ZZZZ, 13 = 


1 




ZZZZ. /} = 


1(1 



1 10 
Time, T 



Figure 5: (Color online) Fidelity error (1 — /) of |0) as a 
function of time for XYXY and ZZ sequences with different 
scales of environmental couplings /3 in units of J. The interval 
is fixed at T = 0.01/ J. Notice that for larger values of /3 the 
data series become almost indistinguishable. The number of 
spin baths ns = 8. 



Clearly, this operation can be implemented by the 
Heisenberg exchange interaction. 

Similarly, the sequence E2 is designed for preserving 
{|EPRi) , IEPR2)} only, and is constructed using Eq. 
(1221) • The cycle for E2 consists of 3 applications of the 
unitary operation (Tp=2,Ds=4 which, upon changing from 
the Bell to to computational basis, takes the form 



/ l+iV3 





E2 



-^+'^ \ 



V 



4 4 
3+i\/3 l-^^/3 
4 4 





-3+1^3 
4 

= exp[!l^(5(i)5(^)-f^a)5p))] 









l-|-iV3 
4 



that is, in terms of an isotropic XZ Hamiltonian. 

Finally, the sequence E3 is designed to preserve all 
four EPR states {|EPRi)}f^g (recaU that it would be 
impossible to preserve exactly 3 orthogonal states in a 
4-dimensional Hilbert space, as the remaining basis el- 
ement would also be inevitably preserved) and is con- 
structed using Eq. (j21[) . Upon transforming, again, to 
the computational basis, the cycle for E3 consists of 4 
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Figure 6: (Color online) Relative survival ^(| V')) (see text) as a 
function of time for various sequences across the Blocli sphere 
equator, with 9 = ±n/2 corresponding to |±X). Notice the 
high survival fidelity near the pointer basis element |0) {0 = 
0). The simulations used ub ~ 8 bath spins with Hb = 0. 

applications of 

i+i 

\^ ^; 

Notice that all the sequences El, E2, and E3 require two- 
body interactions. This is unsurprising since the Bell 
basis is entangled. 

The sequences El, E2, and E3 were numerically simu- 
lated for the two qubits coupled to a bath consisting of 
= 15 spins, and a set of 8 possible initializations (4 
computational basis states and 4 Bell states) was evolved 
for each sequence. The simulation results for El, E2, 
and E3 are summarized in Figs. [3 El and [31 respec- 
tively. Since in this case the contribution to Hsb.c due 
to the internal exchange Hamiltonian Hs need not be 
fully removed by the applied control, we have also plot- 
ted the purity loss (1 — Trp|) to factor out any remaining 
coherent evolution in the system and focus only on de- 
coherence effects when needed. In all these figures, the 
eventual differentiation of the designated PSs from the 
others is clearly visible for each sequence. 

Interestingly, as Fig. [9] reveals, the E3 sequence has 
a side-effect UK ^ 0. We note that the intra-qubit 
Heisenberg coupling Jig has the Bell states as its eigen- 
states. This implies that the E3 sequence effectively 



leaves Hs invariant while it removes the terms in Hsb 
which would otherwise spoil the eigenstates. The net 
effect on the Bell states is as expected, they are pre- 
served as PSs. However, the Heisenberg interaction im- 
plements a nontrivial swap-like evolution in the subspace 
spanned by {|01) , |10)}, while it acts as identity on the 
span of {|00) , |11)}. As a result, the E3 sequence en- 
hances the internal system dynamics by removing the 
unwanted components of Hsn, and effectively resulting 
in a simple "logical gate" [13, ■ This explains the oscil- 
lations in fidelity of {|01) , |10)} that are observed most 
prominently with the E3 sequence. These oscillations are 
absent when we focus on state purity (bottom panel in 
Fig. [5]). Fig. Iini further highlights the difference in the 
evolution of fidelity and purity by explicitly contrasting 
the behavior of computational basis states in the case of 
interacting vs. non-interacting qubits, qualitatively con- 
firming the above picture. 



C. Higher Order Sequences 

The protocols introduced in Section IIVI and quantita- 
tively illustrated in the above examples are all based on 
cancellation of the first order (time-linear) terms in the 
Magnus expansion for the effective cycle Hamiltonian He- 
Higher-order terms in Magnus expansion can also be re- 
moved. Cancellation of higher-order terms in the Magnus 
expansion is wi dely used in perturbative high-order DD 
procedures [3, [3 ■ The timing and pulse type pat- 
terns used in such advanced DD schemes can be adapted 
to the task of PS-engineering if desired. While a de- 
tailed analysis is beyond our current scope, cancellation 
of higher-order Magnus corrections will result, for small 
enough control intervals, in further reducing Hpcr and 
thus the fidelity loss of the preserved PSs. 

More concretely, consider a single qubit, where instead 
of using the ZZ sequences of Sec. IV A[ we apply the 
reflection Q {= (Jz) follows the timing of Uhrig DD [l5| . 
For higher-dimensional systems and for preserving a sin- 
gle PS, the protocol described in Sec. IIV A II can likewise 
be modified to use a control cycle in which (as opposed 
to using two equal length pulse intervals) the pulse inter- 
vals are described by the Uhrig pattern [s^- For more 
than one PS, Uhrig DD is more difficult to adapt to PS- 
engineering. A possible construction may be devised in 
principle for Ds = 2™-dimensional Hilbert spaces, by 
adapting the results in (67j . 

VI. THE ROLE OF IMPERFECT CONTROL 

Throughout the discussion so far, we have assumed 
control resources to be perfect, allowing for precise ini- 
tialization of the system in (one of) the intended PS as 
well as exact implementation of all the required control 
operations. We now revisit these assumptions and ana- 
lyze different ways in which limited control resources can 
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Figure 7: (Color online) Fidelity error (1 — /) and purity loss (1 — Tr(p|)) for the El sequence designated to preserve |EPRi) 
(see text). Different curves (some overlapping) correspond to initialization in the computational basis states and the Bell basis 
states, respectively. The number of bath spins ns = 15, with P — J — K , and the pulse interval t — 0.01/ J. 



impact and/or modify the PS engineering problem. 



A. Preparation Errors 

A deviation of the initial state from the intended 
pointer basis will propagate smoothly to the eventual fi- 
delity of the system. Consider for concreteness a qubit, 
in which the cr^ eigenstates are the designated PSs, and 
let the initial preparation in the PS basis be described by 



1-5A 5{B + iC) 
^0 ~ 1 5{B- iC) 5 A 



(33) 



where B, and C are real numbers and the parameter 
5 quantifies the error strength in the preparation of |0). 
For simplicity, let us further assume that |0) , |1) arc pre- 
served with the maximal fidelity 1 , whereas eigenstates of 
Gx and (J J, evolve to a maximally mixed state as T — cx3. 
Using the linearity of quantum mechanics, we can show 
that at T — > oo, the density matrix is given by 



Poo = 



1-5A 
6A 



The fidelity loss in the evolution of po to po 
leading order in 5 is given by: 

1 



-5\B' + C') + 0{5''). 



up to the 



(34) 



Eq. p4l) implies that a state prepared in a sufficiently 
small neighborhood of a PS (in the convex set of possible 
states, including mixed ones) will evolve with a fidelity 
close to the maximal fidelity. In other words, starting 
with a slightly misprepared initial state and evolving un- 
der a PS-preserving protocol docs result in lower fidelities 
but not a complete fidelity loss. In dynamical systems 
language, the PSs (and all the mixed states diagonal in 
the PS basis) are Lyapunov stable (32l.[33|. The fact that 
the fidelities for states prepared near PSs remains lower 
than the PS-fidelities implies, however, that the dynam- 
ics is not attractive, consistent with the purely unitary 
nature of the applied control. 



B. Pulse Imperfections in Pointer-State Sequences 

While a variety of imperfections can plague control 
Hamiltonians in real experiments, systematic pulse errors 
remain an important limiting factor for the achievable fi- 
delities. As we now show, the latter can be incorporated 
in our theoretical framework as long as they result in a 
constant cycle propagator over the entire control dura- 
tion. Let us reconsider a control cycle Pi , • • • , P„ desig- 
nated for PS-engineering. Let U p. denote the propagator 
associated with the net evolution of the open system dur- 
ing (start to finish) the pulse Pi. The unitary operators 
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Figure 8: (Color online) Fidelity error (1 — /) and purity loss (1 — Tr(p|)) for the E2 sequence designated to preserve 
{|EPRi) , IEPR2)} (see text). Different curves (some overlapping) correspond to initialization in the computational ba-sis 
states and the Bell basis states, respectively. The number of bath spins ub = 15, with P = J = K , and the pulse interval 
T = 0.01/J. The apparently missing short-time fidelity /purity values for {|EPRi) , IEPR2)} are too close to 0.0 to be shown 
on the log scale. 



Up^ are now approximations of the "ideal pulses" Pi: 

Up^ ^ P,cM-^Ep^), (35) 

where Ep. is a Hermitian system-bath operator (the so- 
called "error action" [g^) that we refer to simply as pulse 
error. Following the notation of Sec. IIV A[ instead of Eq. 
([TB| for ideal pulses, the cycle propagator thus reads 

Uc = P„ exp(— ii?p^) exp(— ir„iJ) 
X • • • Pi exp(— lE'pJ exp(— iTiiJ). 

The error model described by Eq. p5p is general enough 
to encompass a large class of imperfections, such as finite- 
width-pulse errors (where Ep. is a system-bath Hamilto- 
nian depending on the pulse implementation) or, even in 
the narrow-pulse limit, rotation-angle and/or rotation- 
axis errors (with Ep. acting on the system only). The 
pulse errors Ep. are systematic in the sense that they only 
depend on the pulse Pi. This results in a fixed (constant) 
cycle propagator over the whole control duration, making 
it still meaningful to use an effective cycle Hamiltonian 
He up to the total time T. Letting s = maxj \\Epj \ \ and 



using the Magnus expansion, we may write [cf. Eq. (fT9|) ] : 

n n 

i=i i=i 
4- 0{T^H^) + 0{nT,Hoe)+0{n^e^), (36) 

where we have used Ej for the error associated with the j- 
th pulse in the sequence, and we have explicitly shown the 
dominant linear terms while denoting the higher-order 
contributions by asymptotic expressions. 

The protocols for engineering PSs discussed in Sec. |IV] 
are all based on bringing He to obey the PS-condition of 
Eq. ([3]), where J^j 'TjQ^HoQj plays the role of ifdom and 
all the remaining terms are gathered into i^per- Clearly, 
the errors associated with pulse imperfections {Ep. ) gen- 
erally belong to Hp^r- Thus, while the task PSl (ensur- 
ing PSs) can be achieved as before, the task PS2 (en- 
suring higher fidelity PSs) will depend on the quality of 
the pulses applied, as well as on access to shorter pulse 
intervals and/or more advanced control cycles. Recall, 
however, that in the protocols described in Sec. IIVI 
a single pulse type {(Jp^Os) is used throughout the se- 
quence, therefore only a single error per pulse appears 
in Eq. Ep^ = Ep for aU j = I,-- - Re- 

markably, for such single-pulse-type protocols, the term 
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Figure 9: (Color online) Fidelity error (1 - /) and purity loss (1 - Tr(p|)) for the E3 sequence designated to preserve the 
complete Bell basis (see text). Different curves (some overlapping) correspond to initialization in the computational basis states 
and the Bell basis states, respectively. The number of bath spins ub ~ 15, with P — J — K , and the pulse interval r = 0.01/ J. 



QjEjQj = Qj^Qj can actually be included in 
Hdom, and is not a perturbative correction. This leaves 

r,Fp„ = T.i/W + 0{nT,Hoe) + 0{n'e'), (37) 

where iJper refers to ideal pulses, and thus results in ad- 
ditional robustness compared to protocols that employ 
different pulse types. 

The crucial implication of the redefinition of i?dom vs. 
H-pci in Eqs. (|36l) and ((37|) is that although this will gen- 
erally imply a worse fidelity lower-bound [cf. Eq. (jlSp ] for 
a PS-protocol implemented with imperfect pulses, such 
imperfections will not result in a degradation of fidelity 
with time. This simple yet practically important point 
is quantitatively illustrated in Fig. [TT] for a single qubit 
interacting with a spin bath as in Sec. IV Al Specifically, 
we have analyzed the combined effect of rotation-angle 
and rotation-axis errors in preserving the state |0) the 
ZZ protocol: the pulse error is characterized by letting 
Ep = r](exCTx + SyO-y + CzOz), where 77 G [0,7r] charac- 
terizes the over-rotation strength, and exTey^ez are ar- 
bitrary (but fixed for all the data points) numbers ran- 
domly sampled from [—1, 1] which characterize the axis 
misalignment. The observed robustness oj the engineered 
PSs with respect to systematic pulse errors is in sharp 
contrast to the behavior of DD protocols where repeat- 
ing control cycles with faulty control tends to consta ntly 
degrade state fidelities for an arbitrarily chosen state [231 ■ 



VII. POINTER STATES FROM IMPERFECT 
DECOUPLING SEQUENCES 

In Sec. IV Al we showed that the imperfect cancellation 
associated with singie-qubit universal DD sequences re- 
sults in the accidental generation of PSs. For instance, 
the XYXY sequence leads to a pointer basis in which the 
Z eigenstates are preserved and, as shown in the previous 
section, can survive without significant fidelity loss even 
in the presence of systematic control errors. Interestingly, 
however, even in situations where decoupling is theoret- 
ically exact, PSs can still emerge from imperfect pulses. 
This manifests in a strongly state-dependent degree of 
stability after many DD cycles, as recently demonstrated 
in the context of DD experiments using the electron spin 
resonance (ESR) of phosphorus donor spins in silicon, see 
in particular the data in Fig. 2 of |22j| and Fig. 1 of [30| . 
In this Section, we show how the emergence of these sta- 
ble states can be naturally re-interpreted and analyzed 
within our general framework, in terms of the emergence 
of a PS-supporting effective Hamiltonian. 

The system under consideration is described in detail 
elsewhere [22|, [s^. For a given spin S', in a frame that 
rotates with the electron's Larmor frequency, the bare 
Hamiltonian describing the interaction between each P 
electron spin and the nuclear spin bath can be approxi- 
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Figure 10: (Color online) Comparison between the evolution of jOl) and |00) under free evolution and the _B3 sequence (see 
text) for K = J (exchange-coupled qubits, Hs 7^ 0) vs. K = (non-interacting qubits, Hs = 0). The number of bath spins 
ub ~ 8, with P = J, and the pulse interval r = 0.01/J. Notice the oscillations of |01) when K ^ 0. 
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Figure 11: (Color online) Fidelity error (1 — /) as a function 
of time for various over-rotation error strengths for the ZZ 
sequence used for preserving |0). The error-free curve {rj = 
0) is practically indistinguishable from the one at low error 
strength {rj = 0.001). The couplings /3 are in units of J, and 
/3 = J. The pulse interval is fixed at r = 0.01/J and the 
number of bath spins ub = 7- 



mated semiclassically approximated by [See Sec. IV A Ij 



where je and bz are the electronic gyromagnetic ratio 
and the total effective magnetic field, respectively. The 
latter accounts for the effect of the nuclear spin environ- 
ment and the off-resonance error resulting from spatial 
inhomogeneity across the sample, both of which can be 
treated as a random variable static in time to a good 
approximation. Clearly, the above interaction Hamilto- 
nian implies a "natural" pointer basis along the Z axis 
in the absence of control: besides the eigenstates of a^, 
other states will not be preserved. In the presence of per- 
fect (instantaneous) pulses, a universal DD cycle (such 
as XYXY or XZXZ) removes Hsb exactly. However, 
in the actual experiments, the DD pulses are not ideal, 
and their imperfections can rapidly accumulate over time, 
eventually nullifying the expected preserving action of 
DD on the spin's quantum state. Despite their destruc- 
tive role in DD, pulse imperfections can in fact select 
PSs, provided that they are constant over time, although 
inhomogeneous across the sample [69| . For ESR experi- 
ments, where the pulses are indeed close to ideal, we can 
approximate them as unitary rotations, with the rota- 
tion axes and angles slightly differing from their nominal 
values. In particular, for X and Y pulses, following the 
notation of Refs. [12, ll^l the actual evolution operators 



Hsb = leSzbz, 
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describing the imperfect rotations are 

Ux = cxp [-i(7r + ej;)(i(T • n)], 
Uy = exp [-i(7r + ej^)(^(T • m)], 

where £x{y)j with ~ £x = are small errors 
in the rotation angles, and n, rfi are the actual rota- 
tion axes slightly differing from their nominal direc- 
tions along X and y, respectively, that is, we may 

write n = { y^l — 77^71^ — v'^n^, V'^y , V^z ) and m = 

{rjiTix, \/l ~ rfrnl. — rfml, mz)^ with 77 ^ 1. The pa- 
rameters e and rj thus play the role of the perturbative 
parameter e in the previous sections. The simplified but 
realistic one-dimensional model derived for ESR on a Si:P 
sample [22|, [s^ suggests that the probability distribu- 
tion for the rotation angle error is peaked around a value 
eo > and is given by P(e) = (l/2eo)[3(l - e/eoT^/^ 
with — 2eo < e < eo- The axis offsets rjny(^z) ^-nd 'r]fTi^(^z) 
each follow a similar distribution peaked around uq > 0. 

Using Eq. ([5]), we can formally define the effective cy- 
cle Hamiltonian He associated with applying the XYXY 
sequence at the inter-pulse interval r (with Tc = 4t). By 
direct calculation, taking into account terms up to the 
second order in pulse errors, we obtain: 

TcHc^ (Jz[-2ri[m.x + riy) + (e^/2) cos fe^r] 

- 2r]{mx + ?iy){[(e/2)(l + sin^^r) - ijiJz cosbzT]ax 
+ [rjniz - (e/2) cosb^T - rjUz smbzT]ay}. 

From the above expression, one can see how the separa- 
tion of He into the dominant and the perturbative term 
happens. As long as off-axis errors are present (77 7^ 0, 
independently of whether rotation errors are present 
also), the dominant term is simply the first contribution, 
TcH^am ~ {rrix + Uy) az, with all remaining terms 
belonging to the perturbative Hamiltonian iJper, which 
is 0[77 max(e, 77)] and is not in any particular direction. 
Using the results of Sec. IIII Al we expect the fidelity loss 
of PSs {Z eigenstates) to be proportional to the square of 
the ratio of the perturbative to dominant Hamiltonian, 
and hence controlled by 0[max(e, 7^)]. In the case where 
77 = 0, to the leading order in the perturbative param- 
eters we have instead TcHc ~ {e^ / 2) (cos bzT)az- Thus, 
while it is formally necessary to redefine i?dom so that it 
includes second-order contributions, the generated PSs 
will still be the Z eigenstates. 

The XZXZ protocol is equivalent to the XYXY one 
in the absence of pulse errors but when pulse errors are 
taken into account, the resulting effective cycle Hamilto- 
nian becomes qualitatively different: 

TcHc^ o-.y[(e/2)(l - sinbzT) - rjnz{l - cos6^r)] 
— [(e/2)(l — sinfe,T) — 77772(1 — cosbzT)] 
X {-7777^2;cr2; + [(e/2)(l + cos 62T) - 77777-2 
-I- 77712 sin 62 t] 0-2}, 



where we took into account that Z pulses in ESR experi- 
ments are implemented as Z = XY ,hy using two closely 
spaced X and Y pulses [22|. Clearly, the dominant term 
is the one in the first row of the above equation, which in- 
cludes only terms of the first order in the pulse errors, and 
implies that the Y eigenstates now emerge as the pointer 
basis. The perturbative correction is 0[max(e, 7y)^] and is 
not in any particular direction. Again, the fidelity loss is 
proportional to the ratio of the perturbative to dominant 
Hamiltonian, and controlled by 0[max(77, (5)]. 

This analysis parallels the consideration of Sec. IV Al 
Thus, it is not surprising that the qualitative conclusions 
about the character of the generated PS are in perfect 
agreement with the experimental and theoretical results 
reported in [2^ [soj . A quantititative analysis would need 
to take into account additional factors, most importantly 
the need to average the fidelity bound over the probabil- 
ity distributions of the offsets and over-rotation errors 
[similar to Eq. ([8])]. While additional technical analy- 
sis of these specific experiments goes beyond our current 
scope, the above clearly demonstrates how our general 
framework can be modified to include the PS generated 
by DD pulse imperfections in context of direct experi- 
mental relevance. 



VIII. CONCLUSION 

We have provided general open-loop unitary control 
protocols for engineering the interaction between a finite- 
dimensional system and its environment in the non- 
Markovian limit, in such a way that a desired set of pure 
states can be maintained as effective pointer states (PSs) 
of the dynamics. Our constructive results are supported 
by analytical upper bounds for the fidelity loss of the en- 
gineered PSs, which allow for systematic improvement by 
simply employing faster and/or more elaborated control 
schemes. While similar in flavor to dynamical decou- 
pling protocols for protecting arbitrary quantum super- 
positions, the methods presented and analyzed in this 
paper aim to provide selective energetic protection for a 
designated set of pure states and their convex combina- 
tions only, effectively allowing for the on-demand gener- 
ation of a robust classical memory. The fact that PSs are 
protected via relative energy gaps in the interaction with 
the environment is manifest in the nature of our perfor- 
mance bounds, which rely on von Neumann mean ergodic 
theorem and simple perturbation theory. Physically, the 
key requirement is to synthesize an effective Hamiltonian 
with a desired "dominant" symmetry structure, and to re- 
duce the "perturbing interactions" that would otherwise 
destabilize the state and thus reduce the fidelity. 

At the expenses of making the control design state- 
dependent, the PS- engineering protocols we have intro- 
duced have distinctive advantages over general-purpose 
dynamical decoupling schemes, which become most 
transparent in the task of engineering complex quantum 
states in multi-qubit systems as PSs. While the task 
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of Bell-state engineering we have analyzed in depth pro- 
vides a paradigmatic example in this respect, the abil- 
ity to engineer arbitrary PSs can potentially be useful 
as a resource in quantum information processing and/or 
quantum metrology tasks. It is also worth stressing that 
the methods we have presented are directly accessible 
via pulsed control, but can modified in principle to al- 
low continuous-time protocols that still effectively trans- 
form the interaction with the environment. Ultimately, 
stronger controls or faster pulse rates are the basic re- 
sources we leverage in achieving high-fidelity PSs within 
our Hamiltonian open quantum system setting. 

A number of further questions may be worth address- 
ing. For instance, it may be interesting to examine 
whether the present framework can be extended to en- 
compass the more general symmetry conditions that al- 
low for time- dependent PSs, as recently considered in [70j . 
From a control-standpoint, the latter problem might in 
turn relate to the possibility of robust time-dependent 
state-tracking, rather than long-time state-preservation 
as examined here. Interestingly, the distance of a density 
operator has recently be invoked to quantify "quantum- 
ness" in the context of coherent energy transfer in bio- 
logical systems [7l| . It is both natural and intriguing to 
ask whether a perturbative mechanism similar to the one 
involved in the generation of stable PSs examined here 
may be brought to bear on the problem of further un- 
derstanding long-lived quantum coherences in complex 
systems. Lastly, the energetic protection against deco- 
herence enjoyed by the PSs may be ultimately connected 
to quantum noiseless subsystem codes. Thus, a natural 
(although possibly highly non-trivial) direction for explo- 
ration is whether the idea of an engineered fidelity guar- 
antee can be viable for a subsystem code that can pre- 
serve genuine quantum information, as opposed to only 
classical information as for a PS basis. 



exists and satisfies [Xll,[7] = [tI]. Note that con- 
vergence is implied in any inner product distance be- 
tween operators that is invariant under the adjoint ac- 
tion X WXU. In the above form, the MET has also 
been invoked to formally relate DD to the quantum Zeno 
effect [ll. 

To visualize the theorem for a finite-dimensional sys- 
tem, consider the complete basis of eigenstates {|</>j;)} of 
U with eigenvalues e"'^' . The operator X can be written 
in the basis of the Eij = operators. The latter 

are transformed according to 



in{(pi—4 



Consequently, as long as there is no degeneracy (4>i ^ 4>j 
unless i = j), we have 



1 1 1 

N ^ ^ N I 

n=0 



-E,,. 



(Al) 



where in the limit of — >■ oo. the r.h.s. approaches zero 
when i j and 1 otherwise. Thus, only diagonal basis el- 
ements are preserved under averaging and the remaining 
ones are annihilated. These diagonal basis elements span 
operators that commute with U. Notice that the left-over 
terms with i ^ j in the limit of a shrinking minimum gap, 
= mini^j — scale with l/{Nfl). 



Appendix B: Initial Fidelity Decay 

While the goal of our scheme is long-time manipulation 
of coherence using periodically repeated control cycles, 
we can nonetheless approximate the initial short time 
behavior of the system. In what follows we assume that 

NS = NTcellH^^Al <C 1, 
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where A^ is the number of control cycles that have been 
applied to the system up to time T = NTc- This allows to 
use the linear (in S) component of the interaction picture 
propagator exp(— iJlJ^' ) as a substitute for the propagator 
when PS fidelity is concerned. We can estimate the norm 
of the effective Hamiltonian H^' even before using Eq. 
(|Aip . directly from Eq. pT|) in the main text: 



\n^;,'\\<N\\E\\ 



NWHr, 



N6. 



Appendix A: Mean Ergodic Theorem 

Let U and X denote, respectively, a unitary opera- 
tor and arbitrary bounded operator acting on the same 
Hilbert space. The von Neumann's Mean Ergodie Theo- 
rem (MET) implies that the limit 



lim — 

N^ao N 



N-1 



J2 u-'^xu"" = a" 



n=0 



The fidelity loss of the system after N cycles can 
be bounded starting from the following bound for the 
Uhlmann fidelity = fj^^: 



1- <d[\0){0\, exp(-zr!W ) |0) (0| expi+iQ^^^) 
where D[pi,p2\ is the trace distance: 

D[pi,P2] = ^IIpi - P2|li. 
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We can then use (1 — /at) < 2(1 — /^) [following from 
l — x< 2(1 — x^/^) for < a; < 1] and the general bounds 
in Ref. [S^ , to show that 

1 - In < cxpiwn^;:,^ + 0{nH^)\\) - I 

<(e-l)||l]W||<(e,-l)iV5, (Bl) 

where in the last step we used — 1 < (e — l)a;, for 
a; < 1. We have thus bounded the short-time fidelity 
loss in terms of a linear function of the number of cycles 
N for small N. Note that for a single cycle, a tighter 
fidelity bound can be given, directly in terms of the (bath- 
averaged) variance of -ffper in the initial PS [1^, leading 
to 1 - /(Tc) < (e||7Jper||?T2. For large values of N, the 
results of Sec. |in]are applicable instead. 

Appendix C: Implementing Arbitrary Reflections 

Implementing the PS-preserving pulses defined in Sec. 
IIVI is straightforward provided that tunable control 
Hamiltonians diagonal in the pointer basis are available, 
that is, Hamiltonians of the form 

Ds-l 
i=0 



Note that in principle, even a constant in time Hamil- 
tonian which is diagonal in the pointer basis and has 
non-commensurate eigenvalues would suffice for generat- 
ing any unitary which is diagonal in the pointer basis [tJ . 
If, however, no such control Hamiltonians are available, 
the problem of producing the required pulse operators 
can be difficult. 

In particular, the pulse operators used in our proto- 
cols can be entangling if the system is multipartite. For 
example, this was explicitly the case in the PS-protocols 
employed for Bell-state engineering. Sec. IV Bl Further- 
more, note that the operator Q used in Sec. IIV A II to 
preserve a single PS corresponds to a multiple-controUed- 
phase quantum gate, thus, generally, entangling opera- 
tions are required also to preserve a single product state 
(say, an element of the computational basis). Interest- 
ingly, however, preserving all computational basis ele- 
ments in a multiple-qubit system requires no entangle- 
ment. For example, consider the control cycle resulting 
from concatenating (nesting) the ZZ cycles acting on ev- 
ery qubit jisl \q1\ , resulting in preservation of the compu- 
tational basis as an engineered pointer basis. These ideas 
suggest possible connections to the quantum search algo- 
rithm [75| , in which the oracle performs the same (entan- 
gling) operation Q to select a particular element in the 
computational basis. 
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